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1. Warm-up: basic facts about s dvv and b ^ sy 

The s ^ dvv transition is one of the rare examples of an electroweak (EW) process whose 
leading contribution starts at 0{Gp) within the standard model (SM). At the one-loop level it 
proceeds through Z-penguin and EW box diagrams which are highly sensitive to the underlying 
short-distance (SD) dynamics. Sample diagrams are shown on the left of Fig. 1. Separating the 
contributions according to the intermediate up-type quark running inside the loops, the QCD cor- 
rected amplitude takes the form 

Asm(^-^vv)= £ v;v^,x|^oc^(A5 + a5) + ^ln^A + ^A, (1.1) 

q=u,c,t w "'-VI "'vi ■'"h' 

where Vij denote the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix and A = | V„i| = 
0.225. The hierarchy of the CKM elements would obviously favor the charm and up quark con- 
tributions, but the power-like Glashow-Iliopoulos-Maiani (GIM) mechanism, arising mainly from 
the SU {2)l breaking in the Z-penguin amplitude, leads to a very different picture. The top quark 
contribution, caiiying a large CP-violating phase, accounts for 68% of Asm(^ dvv), while cor- 
rections due to internal charm and up quarks amount to ~ 29% and a mere ~ 3%. These properties 
imply that long-distance (LD) effects in the direct CP-violating Ki n^vv mode ai^e negligible, 
while they are highly suppressed in — > 7r+vv(7). 

A related important feature, following from the EW structure of AsmC'^ —>■ dvv) as well, is that 
the SD contributions to both K ^ nvv decay modes are governed by a single effective operator, 
namely 

Qv = {sLYndL){vLY^v,). (1.2) 

By virtue of the conservation of the V — A current, large QCD logarithms appear only in the charm 
quark contribution to Eq. (1.1). The hadronic matrix element of Qy itself, can be extracted very 
precisely from the wealth of available data on ^ — > n£v (Ki^^) decays. 

In summary, the superb theoretical cleanness and the enhanced sensitivity to both non-standard 
flavor and CP violation, make the s ^ dvv channels unique tools to discover or, if no deviation 
is found, to set severe constraints on non-minimal-flavor-violating (MFV) physics where the hard 
GIM cancellation present in the SM and MFV is in general no longer active [1]. 

Unlike s — > dvv, the b ^ sy transition is dominated by perturbative QCD effects which re- 
place the power-like GIM mechanism present in the EW vertex by a logarithmic one. The mild 
suppression of the QCD corrected SM amplitude 

AMb^sy)= £ v;,V,.Kl^^\n'^X^ + \A^+\n'^X\ (1.3) 

q=u.c.t '-^'^ ^"IV 

reduces the sensitivity of the process to high scale physics, but enhances the B — > XsY branching 
ratio (BR) with respect to the purely EW prediction by a factor of around three. The logarithmic 
GIM cancellation originates from the non-conservation of the effective tensor operator 

Qi = i^mh{hOy,b,)F^\ (1.4) 

which is generated at the EW scale by photon penguin diagrams involving W-boson and top quark 
exchange. Sample one- and two-loop diagrams ai^e shown on the right of Fig. 1. 
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Figure 1: Feynman diagrams contributing to the s dvv (left) and b — > sy (right) transition in the SM. 
Their scaling behavior in the heavy top quark mass limit is also shown. See text for details. 



After including logarithmic enhanced QCD effects, the dominant contribution to the partonic 
b ^partonic^ dccay rate stems from charm quark loops that amount to ~ 158% of Asuib sy)- 
The top contribution is with ~ —60% of Asm(^ —>■ ^7) compared to the charm quark effects less 
than half as big and has the opposite sign. Diagrams involving up quarks are suppressed by small 
CKM factors and lead to an effect of only ~ 2% in Asm(^ —>■ ^7)- Due to the inclusive character 
of the B XsY mode and the heaviness of the bottom quark, » A ~ Aqcd> non-perturbative 
effects arise only as small corrections to the partonic decay rate. 

In summary, the good theoretical control on and experimental accessibility of the b ^ sy rate, 
and its large generic sensitivity to non-standard sources of flavor and CP violation, allows to derive 
stringent constraints on a variety of new physics (NP) models, in particular on those where the 
flavor-violating chiral transition of the amplitude is not suppressed. 



2. Recent theoretical progress in iS^^ ^ :;r^vv and ^ 'JI^vv{y) 

After summation over the three lepton families the SM BRs of the K — > nvv modes read 

^r{Kl^k\v)sm = Kl(^^X^ , (2.1) 
BR(^+ ^ 7r+vv(7))sM = '^+(1 + Arm) 



, (2.2) 



where A,- = V*Vid. The top quark contribution X = 1 .456 ± 0.017^, ± 0.013;;, ± 0.015ew is known 
through next-to-leading order (NLO) in QCD [2]. Its overall uncertainty of ~ 2% is in equal shares 
due to the parametric error on the top quark mass, the sensitivity on the matching scale fit, and 
two-loop EW effects for which only the leading term in the heavy top quark mass expansion has 
been calculated [3] . 

Major theoretical progress has been recently made concerning the extraction of the hadronic 
{7l'\sY^d\Kj) matrix elements from Kci, data, by extending the classic chiral perturbation the- 
ory (ChPT) analysis of leading order (LO) 0{p^e^^^) isospin-breaking effects [4]. Here e^^' oc 
(m„ — md)/ms. The inclusion of NLO 0{p^e^^^) and partial next-to-next-to-leading order (NNLO) 
0{p^e^^^) corrections in the ChPT expansion [5] leads to a reduction of the uncertainties on the 
Kl Ti^vv and —>■ tt+vv matrix elements by a factor of 4 and ~ 7. Since the overall uncer- 
tainties on Kl = (2.229 ±0.017) x lO^^^ (A /0.225)** and K+ = (5.168 ±0.025) x lO^i" (A /0.225)^ 
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Figure 2: Charm quark mass renormalization scale jic dependence of Pc (left) and BR(B Xs7)sm (right) 
at LO, NLO, and NNLO in QCD. The width of the curves on the left indicate the uncertainty due to higher 
order terms in that affect the evaluation of (Xs{lJ.c) from as{Mz)- See text for details. 

[5] are now dominated by experimental errors, a further improvement in the extraction of the rare 
.S'-decay matrix elements will be possible with improved data for the K/t, slopes and Kj^ BRs. LD 
QED corrections affecting 7r+vv(7) are encoded by Aem in Eq. (2.2). At LO in the ChPT 

expansion, these infrared finite 0{p^a) corrections amount to Aem = —0.003 [5] for a maximum 
energy of 20MeV of the undetected photon. More details on the ChPT analysis of rare /T-decay 
matrix elements can be found in [5] . 

The parameter Pc entering Eq. (2.2) results from Z-penguin and EW box diagrams involv- 
ing internal charm quark exchange. As now both high- and low-energy scales are involved, a 
complete renormalization group analysis of this term is required. In this manner, large loga- 
rithms \nmc/M„ are resummed to all orders in a^. The inclusion of NNLO QCD corrections 
[6] leads to a significant reduction of the theoretical uncertainty by a factor of ~ 4, as it removes 
almost the entire sensitivities of Pc on the charm quark mass renormalization scale pLc and on 
higher order terms in that affect the evaluation of as{pic) from av(Mz). This is illustrated by 
the plot in the left panel of Fig. 2. For nic = (1.30 ± 0.05) GeV one obtains the NNLO value 
Pc = (0.374 ± 0.03 1,„^ ±0.009pert ±0.009aJ (0.225/A)4 [6], where the individual errors are due to 
the uncertainty on the charm quark MS mass, higher-order perturbative effects, and the parametric 
error on as{Mz) . Since the residual error on Pc is now fully dominated by the parametric uncertainty 
coming from nic, a better determination of the charm quark mass is clearly an important theoretical 
goal in connection with — > 7r+vv(7). 

To gain an accuracy on the 7r+vv(7) BR of a few percent, it is necessary to account for 

subleading effects not described by the effective Hamiltonian that includes the dimension-six oper- 
ator Qy of Eq. (1.2). The subleading corrections can be divided into two groups: /) contributions of 
dimension-eight four fermion operators generated at the charm quark mass renormalization scale 
pic [7,8], and //) genuine LD contributions due to up quark loops which can be described within the 
framework of ChPT [8]. Both contributions can be effectively included by SPc^u = 0.04 ±0.02 [8] 
in Eq. (2.2). Numerically, they lead to an enhancement of BR(A'+ 7r+vv(7))sM by ~ 7%. The 
quoted residual error of 5Pc,u could in principle be reduced by means of a dedicated lattice QCD 
computation [9]. 

Taking into account all the indirect constraints from the latest unitarity triangle fit [10], one 
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Figure 3: Error budget of the SM prediction of BR(ii:i n^vv) (left) and BR(7i:+ n+vv{Y)) (right). 
See text for details. 

finds, by adding errors quadratically, the following SM predictions for the two K —>■ %vv rates 

V,^{Kl 7r°vv)sM = (2.54±0.35) x 10"" , (2.3) 
BR(i^+ ^ 7r+vv(7))sM = (7.96±0.86) x 10"" . (2.4) 

The error budgets of the SM predictions of both ^ — > TTVV decays are illustrated by the pie charts 
in Fig. 3. As the breakdown of the residual uncertainties shows, both decay modes are at present 
subject mainly to parametric errors (81%, 69%) stemming from the CKM parameters, the quark 
masses nic and m,, and a^^Mz). The non-parametric errors in Ki^ —>■ n^vv are dominated by the 
uncertainty due to higher-order perturbative effects (15%), while in the case of —>■ 7r+vv(7) the 
errors due to dimension-eight charm and LD up quark effects (16%) and left over scale uncertainties 
(12%) are similar in size. Given the expected improvement in the extraction of the CKM elements 
through the B-factories, SM predictions for both K Tivv rates with an accuracy significantly 
below 10% should be possible before the end of this decade. Such precisions are unique in the field 
of flavor-changing-neutral-current processes. 

3. Recent theoretical progress in 5 ^ XsY and B Xsi^i~ 

Considerable effort has gone into the calculation of fixed-order logarithmic enhanced NNLO 
QCD corrections to B ^ X^y [11, 12, 13, 14]. A crucial part of the NNLO calculation is the 
interpolation in the charm quark mass performed in [13]. The three-loop 0{aj) matrix elements 
of the current-current operators Q12 contain the charm quark, and the NNLO calculation of these 
matrix elements is essential to reduce the overall theoretical uncertainty of the SM calculation. In 
fact, the largest part of the theoretical uncertainty in the NLO analysis of the BR is related to the 
definition of the mass of the charm quark [15] that enters the 0{as) matrix elements {sY\Qi2\b). 
The latter matrix elements are non- vanishing at two loops only and the scale at which trie should be 
normalized is therefore undetermined at NLO. Since varying nic between mc{mc) ~ 1.25GeV and 
mc{mb) ^ 0.85 GeV leads to a shift in the NLO BR of more than 10% this issue is not an academic 
one. 

Finding the complete NNLO correction to {sy\Q\2\b) is a formidable task, since it involves 
the evaluation of hundreds of three-loop on-shell vertex diagrams that are presently not even known 
in the case = 0. The approximation made in [13] is based on the observation that at the phys- 
ical point rric ~ 0.25 Mb the large m^. S> nib asymptotic form of the exact 0{as) [16] and large-jSo 
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0{ajpo) [12] result matches the small ^ nih expansion rather well. This feature prompted the 
analytic calculation of the leading term in the » mh expansion of the three-loop diagrams, and 
to use the obtained information to perform an interpolation to smaller values of assuming the 
0{ajPo) part to be a good approximation of the full 0{aj) result for vanishing charm quark mass. 
The uncertainty related to this procedure has been assessed in [13] by employing three ansatze 
with different boundary conditions at = 0. A complete calculation of the 0{aj) corrections to 
(■^7121,21^) in the latter Umit or, if possible, for nic ~ 0.25 nifo, would resolve this ambiguity and 
should therefore be attempted. 

Combining the aforementioned results it was possible to obtain the first theoretical estimate of 
the total BR of B ^ X^y at NNLO. For a photon energy cut of Ey > Ecut with E^ut = 1 .6GeV in the 
B-meson rest-frame the improved SM evaluation is given by [13, 14] 

BR(B^X,7)sM = (3.15 ±0.23) x 10"^ (3.1) 

where the uncertainties from hadronic power corrections (5%), parametric dependences (3%), 
higher-order perturbative effects (3%), and the interpolation in the charm quark mass (3%) have 
been added in quadrature to obtain the total eiTor. 

The reduction of the renormalization scale dependence at NNLO is clearly seen in the right 
panel of Fig. 2. The most pronounced effect occurs in the case of He that was the main source of 
uncertainty at NLO. The cun^ent uncertainty of 3% due to higher-order effects is estimated from 
the variations of the NNLO estimate under change of renormalization scales. The central value in 
Eq. (3.1) corresponds to the choice IJ.w,b,c = (160,2.5, 1.5) GeV. More details on the phenomeno- 
logical analysis including the list of input parameters can be found in [13]. 

It is well-known that the operator product expansion (OPE) for B —>■ X^y has certain limita- 
tions which stem from the fact that the photon has a partonic substructure. In particular, the local 
expansion does not apply to contributions from operators other than Qj, in which the photon cou- 
ples to light quarks [17, 18]. While the presence of non-local power corrections was thus foreseen 
such terms have been studied until recently only in the case of the (28,28) interference [17]. In 
[19] the analysis of non-perturbative effects that go beyond the local OPE have been extended 
to the enhanced non-local terms emerging from {Qi,Q^) insertions. The found correction scales 
like 0{asA/mh) and its effect on the BR was estimated using the vacuum insertion approximation 
to be —[0.3,3.0]%. A measurement of the flavor asymmetry between — > Xgy and B Xsy 
could help to sustain this numerical estimate [19]. Potentially as important as the latter corrections 
are those arising from the (21.2,2?) interference. Naive dimensional analysis suggests that some 
non-perturbative corrections to them also scale like 0{asA/mi,). Since at the moment there is not 
even an estimate of those corrections, a non-perturbative uncertainty of 5% has been assigned to 
the result in Eq. (3.1). This eiTor is the dominant theoretical uncertainty at present and thought 
to include all known [19] and unknown 0{asA/mh) terms. Calculating the precise impact of the 
enhanced non-local power corrections may remain notoriously difficult given the limited control 
over non-perturbative effects on the light cone. 

A further complication in the calculation of B ^ X^y arises from the fact that all measurements 
impose stringent cuts on the photon energy to suppress the background from other B-meson decay 
processes. Restricting Ey to be close to the physical endpoint Smax = 'Mb/2, leads to a breakdown 
of the local OPE, which can be cured by resummation of an infinite set of leading-twist terms 
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into a non-perturbative shape function [20]. A detailed knowledge of the shape function and other 
subleading effects is required to extrapolate the measurements to a region where the conventional 
OPE can be trusted. 

The transition from the shape function to the OPE region can be described by a multi-scale 
OPE (MSOPE) [21]. In addition to the hard scale jU/, ~ ~ 5GeV, this expansion involves a 
hard-coUinear scale /i/,c V mi, A ~ 2.5 GeV coiTcsponding to the typical hadronic invariant mass 
of the final state Xs, and a soft scale /i^ A ~ 1.5GeV related to the width A/2 = mb/2 — Ecut 
of the energy window in which the photon spectrum is measured. In the MSOPE framework, the 
perturbative tail of the spectrum receives calculable corrections at all three scales, and may be 
subject to large perturbative corrections due to the presence of terms proportional to av(\/'MfoA) 
0.27 and a, (A) ~0.36. 

A systematic MSOPE analysis of the {Qi,Qi) interference at NNLO has been performed 
in [22]. Besides the hard matching corrections, it involves the two-loop logarithmic and con- 
stant terms of the jet [21, 24] and soft function [25]. The three-loop anomalous dimension of 
the shape function remains unknown and is not included. The MSOPE result can be combined 
with the fixed-order prediction by computing the fraction of events \ — T that lies in the range 
£^cut = [1-0, 1.6] GeV. The analysis [22] yields 

1 _ r = 0.07+0;03p^^ ± 0.02hadr ± 0.02pa,s , (3.2) 

where the individual errors are perturbative, hadronic, and parametric. The quoted value is almost 
twice as large as the NNLO estimate \ — T = 0.04ib0.01pert obtained in fixed-order perturbation 
theory [13, 14, 23] and plagued by a significant additional theoretical error related to low-scale 
perturbative corrections. These large residual scale uncertainties indicate a slow convergence of 
the MSOPE series expansion in the tail region of the photon energy spectrum. Given that A is 
always larger than 1 .4 GeV and thus fully in the perturbative regime this feature is unexpected. 

Additional theoretical information on the shape of the photon energy spectrum can be obtained 
from the universality of soft and collinear gluon radiation. Such an approach can be used to pre- 
dict large logarithms of the form ln(£'max — £^cut)- These computations have also achieved NNLO 
accuracy [26] and incorporate Sudakov and renormalon resummation via dressed gluon exponenti- 
ation (DGE) [26, 27]. The present NNLO estimate of 1 - T = 0.016 ± 0.003pert [26, 28] indicates 
a much thinner tail of the photon energy spectrum and a considerable smaller perturbative uncer- 
tainty than reported in [22] . The DGE analysis thus supports the view that the integrated photon 
energy spectrum below = 1.6 GeV is well approximated by a fixed-order perturbative calcula- 
tion, complemented by local OPE power corrections. To understand how precisely the tail of the 
photon energy spectrum can be calculated requires nevertheless further theoretical investigations. 

The study of ft — > sl'^i'^ transitions can yield useful complementary information, when con- 
fronted with the less rare b ^ sj decays, in testing the flavor sector of the SM. In particular, a 
precise measurement of the inclusive B X^i^i^ decay distributions would be welcome in view 
of NP searches, because they are amenable to clean theoretical descriptions for dilepton invariant 
masses in the ranges IGeV^ <q^ < GGeV^ [29, 30, 31] and > UGeV^ [32]. 

The SM calculations of the differential rate and forward-backward (FB) asymmetry have both 
reached NNLO precision [29, 30, 33]. In the case of the the dilepton invariant mass spectrum. 
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integrated over the low-^^ region, the most recent SM prediction reads [31] 

BR(B^X,i+r)^^'^'<'^'<^'^'^' = (1.59±0.11) X 10-^ (3.3) 

The position of the zero of the FB asymmetry is known to be especially sensitive to NP effects. In 
the SM one finds [30] 

^o,SM = (3.V6±0.33)GeV2. (3.4) 

The total errors in Eqs. (3.3) and (3.4) have been obtained by adding the individual parametric and 
perturbative uncertainties in quadrature. Besides the differential rate and FB asymmetry, an angular 
decomposition of S ^ X^i'^l^ provides a third observable that is sensitive to a different combina- 
tion of Wilson coefficients. This recent observation [36] might allow to extract SD information 
from limited data on B ^ X^i^t^ in the lov/-q^ region more efficiently. 

Like in the case of B ^ X^J, experimental cuts complicate the theoretical description of 
B Xst^t~ as they make the measured decay distributions sensitive to the non-perturbative shape 
function [34]. In paiticulai^, putting an upper cut M|^' = [1.8,2.0] GeV on the hadronic invariant 
mass of Xs [35], in order to suppress the background from B Xclv Xst^l'vv at small q^, 
causes a reduction of the rate by (10 — 30)% [34]. Although the reduction can be accurately calcu- 
lated using the universal B — > XsY shape function [34], subleading shape functions may introduce 
an additional error of 5% in Eq. (3.3). Similarly, no additional uncertainty for unknown subleading 
non-perturbative corrections has been included in Eqs. (3.3) and (3.4). Most importantly, uncalcu- 
lated 0{asA/mi,) non-perturbative corrections may imply an additional uncertainty of 5% in the 
above formulas. This issue deserves an dedicated study. 

4. CMFV: combining K+ n+vvir), B XsY, B X^i+r, and Z ^ bb 

A way to study possible NP effects in flavor physics consists in constraining the Wilson co- 
efficients of the operators in the low-energy effective theory. This model-independent approach to 
MFV has been applied combining various K- and B-meson decay modes both including [37] and 
neglecting [38, 39, 40] operators that do not contribute in the SM, that is, so-called constrained 
MFV (CMFV) [41] scenarios. 

The main goal of the recent CMFV study [40] was an improved determination of the range 
allowed for the NP contribution AC = C — Csm to the universal Z-penguin by a thorough global fit 
to the Z ^ bb pseudo observables (POs) 7?^, Aj„ and Apg [42] and the measured B — > XgJ [43] and 
B — > Xsi^i^ [35] rates. The derived model-independent bounds 

AC = -0.026 ±0.264 (68% CL), AC = [-0.483,0.368] (95% CL) , (4.1) 

imply that large negative contributions that would reverse the sign of the SM Z-penguin amplitude 
are highly disfavored in CMFV scenarios due to the strong constraint from [40]. Interestingly, 
such a conclusion cannot be drawn by considering only flavor constraints [39], since a combination 
of BR(B X,7), BR(B X,i+r), and BR{K+ k+vv{y)) does not allow to distinguish the 
SM solution AC = from the wrong-sign case AC ~ —2 at present. The constraint on AC within 
CMFV following from the simultaneous use of R^, Ab, Apg, BR(B X.y), and BR(B XJ+i-) 
can be seen in the left panel of Fig. 4. 
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Figure 4: Left: Constraints on ACf^ and AC within CMFV that follow from a combination of the mea- 
surements of B ^ XsY, B Xst^l^ , and Z bb. The colors encode the frequentist 1 — CL level and the 
corresponding 68% CL and 95%CL regions as indicated by the bars on the right side of the panel. Right: 
Frequentist 1 - CL level of BR{K+ 7t+Vv{Y)) vs. BR{Kl 7t°vv) in CMFV. The two black ellipses 
indicate the 68% CL and 95% CL regions in the SM while the best fit value within CMFV is specified by the 
black cross. The central value of BR(/r+ ^^vv{y))sm as a function of the charm quark MS mass in GeV 
and the lower experimental 90% CL bound and central value are also shown. See text for details. 



One can also infer from this figure that two regions, resembling the two possible signs of the 
amplitude A{b — > sy) Cj^{mh), satisfy all existing experimental bounds. The best fit value for 
the NP contribution ACf = Cf{mb) - Cff^ (mh) is very close to the SM point residing in the 
origin, while the wrong-sign solution located on the right is highly disfavored, as it corresponds 
to a BR(S — > Xsl'^i^) value considerably higher than the measurements [44]. The corresponding 
limits are [40] 

ACf = -0.039 ±0.043 (68% CL) , ACf = [-0.104,0.026] U [0.890,0.968] (95% CL) . (4.2) 

Similar bounds have been presented previously in [39]. Since BR(S Xv7)sm as given in Eq. (3. 1) 
is lower than the experimental world average BR(B X^y) = (3.55 ±0.26) x 10""* [45] by 1.2 a, 
extensions of the SM that predict a suppression of the b ^ sy amplitude are strongly constrained. 
In particular, even the SM point ACy^ = is almost disfavored at 68% CL by the global fit. 

The stringent bound on the NP contribution AC given in Eq. (4.1) translates into tight two-sided 
hmits for the BRs of all Z-penguin dominated K- and B-decays [40]. In the case of the K nvv 
modes the allowed ranges for the CMFV BRs read [40] 

BR(^L ^ 7r°vv)cMFV = [1.55,4.38] x 10"'^ (95%CL) , 
m.{K+ 7r+vv(7))cMFV = [4.29, 10.72] x 10"^^ (95%CL) , ^^'^"^ 

when NP contributions to EW boxes are neglected. These bounds imply that in CMFV models 
the Ki TT^vv and 7r+vv(7) rates can differ by at most and +3q^" from their SM 

expectations. In particular, the Grossman-Nir bound BR(iCx 7r°vv) < 4.4BR(/(r+ k^vv{y)) 
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[46] following from isospin symmetry, cannot be saturated in CMFV models. Tiiis is illustrated 
by the plot in the right panel of Fig. 4 which shows the frequentist 1 — CL level of the prediction 
BR(^+ ^ 7r+vv(7))cMFv v^. BR{Kl 7r°vv)cMFv- 

A strong violation of the CMFV bounds in Eq. (4.3) by future precision measurements of 
BR{Ki K^vv) and/or BR(^+ 7r+vv(7)) will imply a failure of the CMFV assumption, sig- 
naling either the presence of new effective operators and/or new flavor and CP violation. A way 
to evade the above limits is the presence of sizable EW box contributions. While these possibility 
cannot be fully excluded, general arguments and explicit calculations indicate that it is difficult to 
realize in the CMFV framework. 
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